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Abstract 

A review on the application of Melnikov 's method to control homo- 
clinic and heteroclinic chaos in low-dimensional, non- autonomous and dis- 
sipative, oscillator systems by weak harmonic excitations is presented, in- 
cluding diverse applications such as chaotic escape from a potential well, 
chaotic solitons in Frenkel-Kontorova chains, and chaotic charged particles 
in the field of an electrostatic wave packet. 
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1 Introduction 

During the past 15 years or so, diverse techniques of non-feedback chaos control 
have been proposed (Chen & Dong 1998) that may be roughly classified into 
three types: (i) The parametric excitation of an experimentally adjustable pa- 
rameter; (ii) Entrainment to the target dynamics; and (iii) The application of 
a coordinate-independent (or dipole) external periodic excitation. It is shown 
below that techniques (i) and (iii) may be unified in a general setting for the 
class of dissipative systems considered in this present review. There exists nu- 
merical, experimental, and theoretical evidence that the period of the most 
effective chaos-controlling excitations usually is a rational fraction of a certain 
period associated with the uncontrolled system, although the effectiveness of in- 
commensurate excitations has also been demonstrated in some particular cases 
(Chacon & Martmez 2002). Indeed, resonances between the chaos-controlling 
excitation and (i) a periodic chaos-inducing excitation, (ii) an unstable periodic 
orbit embedded in the chaotic attractor, (iii) a natural period in a period dou- 
bling route to chaos, or (iv) a period associated with some peak in the power 
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spectrum, have been considered in diverse successful chaos-controlling excita- 
tions. This is not really surprising since these types of resonances are closely 
related to each other. For instance, when a damped, harmonically forced os- 
cillator exhibits a steady chaotic state, the power spectrum corresponding to 
a given state variable typically presents its main peaks at frequencies which 
are rational fractions of the chaos-inducing frequency for certain ranges of the 
chaos-inducing amplitude. The extensive literature concerning experimental, 
theoretical, and numerical studies of non-feedback methods is frankly unap- 
proachable because of its volume in a review of the present type. Therefore, 
only pioneering key work (from the author's viewpoint) is mentioned in the fol- 
lowing. The effectiveness of periodic parametric excitations in suppressing chaos 
was shown by Alekseev & Loskutov (1987). Hiibler & Liischer (1989) discussed 
how a nonlinear oscillator can be driven toward a given target dynamics by 
means of resonant excitations. Braiman & Goldhirsch (1991) provided numeri- 
cal evidence to show the possibility of inhibiting chaos by an additional periodic 
coordinate-independent excitation. Salerno (1991) showed, by the analysis of a 
phase-locked map, the possibility of suppressing chaos in long biharmonically 
driven Josephson junctions. Chacon & Diaz Bejarano (1993) discussed a new 
way to reduce or suppress steady chaotic states, by only altering the geometrical 
shape of weak periodic perturbations. Kivshar (1994) showed analytically 

and numerically that the suppression of chaos may be effectively achieved by 
applying a high-frequency parametric force to a chaotic dynamical system. Ex- 
perimental control of chaos by weak periodic excitations has been demonstrated 
in many diverse systems, including magnetoelastic systems (Ditto et al. 1990), 
ferromagnetic systems (Azevcdo & Rezcnde 1991), electronic systems (Hunt 
1991), laser systems (Roy et al. 1992; Meucci et al. 1994; Chizhevsky & Cor- 
balan 1996; Uchida et al. 1998), chemical reactions (Petrov et al. 1993; Alonso 
et al. 2003), neurological systems (Schiff et al. 1994), and plasma systems (Ding 
et al. 1994). 

This paper summarizes some main results concerning the application of Mel- 

nikov's method (Melnikov 1963; Arnold 1964; Guckcnheimcr & Holmes 1983; 
Wiggings 1990) to the problem of control of chaos in low-dimensional, non- 
autonomous and dissipative, oscillator systems by small-amplitude harmonic 
perturbations. Specifically, the class of systems considered is described by the 
differential equation 

^+-^ = -d{x,x) +Pc{x,x)Fc{t) +ps {x,x) Fs{t), (1.1) 

where U{x) is a nonlinear potential, —d{x,x) is a generic dissipative force 
which may include constant forces and time-delay terms, Pc{x, x)Fi,(t) is a 
chaos- inducing excitation, and pg (x, i) Fg{t) is an as yet undetermined suitable 
chaos-controlling excitation, with Ff.{t), Fg{t) being harmonic functions of ini- 
tial phases 0, O, and frequencies cj, Q, respectively. It is worth mentioning that 
Melnikov 's method imposes on (1.1) some additional limitations: the excitation, 
time-delay, and dissipation terms are weak perturbations of the underlying con- 
servative system x + dU{x)/dx = which has a separatrix. The original work 
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of Melnikov (1963) was generalized by Arnold (1964) to a particular instance of 
a time-periodic Hamiltonian perturbation of a two-degree-of-freedom integrable 
Hamiltonian system. Fifteen years later, Holmes (1979) was the first to apply 
Melnikov 's method (to a damped forced two- well Duffing oscillator) in the west. 
From then on the method began to be popular. Chow et al. (1980) rediscovered 
Melnikov 's results using alternative methods and emphasized that homoclinic 
and subharmonic bifurcations are closely related. Through the 1980s a great 
variety of extensions and generalizations of Melnikov 's approach were devel- 
oped (Greenspan 1981; Holmes & Marsden 1982; Lerman & Umanski 1984; 
Grcundlcr 1985; Salam 1987; Schecter 1987; Wiggings 1987). The interested 
reader is referred to the books by Lichtenberg & Lieberman (1983), Gucken- 
heimer & Holmes (1983), Wiggings (1988), and Arrowsmith & Place (1990) for 
more details and references. The work of Cai et al. (2002) provides the simplest 
extension of Melnikov 's method to include perturbational time-delay terms. 

The application of Melnikov 's method to controlling chaos in low-dimensional 
systems by weak periodic perturbations began in about 1990. Indeed, Lima & 
Pettini (1990) provided a heuristic argument to extend the idea that parametric 
perturbations can modify the stability of hyperbolic or elliptic fixed points, in 
the phase space of linear systems, to the case of nonlinear systems, and hence 
that parametric perturbations could also provide a means to reduce or suppress 
chaos in nonlinear systems. They used for the first time the Melnikov method to 
analytically demonstrate this conjecture in the case of a damped driven two- well 
Duffing oscillator subjected to a chaos-suppressing parametric excitation. How- 
ever, their insufBcient analysis of the corresponding Melnikov function led them 
into gross errors in their final results and conclusions. Specifically, they failed 
both theoretically to demonstrate the sensitivity of the suppression scenario to 
the initial phase of the chaos-suppressing excitation and to find it numerically. 
They also failed theoretically to predict the suppression of chaos in the case of 
subharmonic resonances (between the chaos-inducing and chaos-suppressing ex- 
citations) higher than the main one. Although a part of their erroneous analysis 
of the Melnikov function originated from a mistake in its calculation (Cuadros 
& Chacon 1993; Lima & Pettini 1993), its main weakness was in not providing 
a correct necessary and sufficient condition for the Melnikov function to always 
have the same sign (i.e., for the frustration of homoclinic bifurcations). For 
the two-well Duffing oscillator that they considered, such a correct necessary 
and sufficient condition was first deduced for the general case of subharmonic 
resonances by Chacon (1995a), where the extremely important role of the ini- 
tial phase (of the chaos-suppressing excitation) on the suppression scenario was 
demonstrated theoretically. Cicogna & Fronzoni (1990) studied the suppression 
of chaos in the Josephson-junction model 

-I- [1 -I- ^ cos {Q.t + 6)] sin = -dij) + 7 cos {ujt) , (1.2) 

where the parametric excitation ^ cos (ili -I- 9) sin (p is the chaos-suppressing ex- 
citation, for the single case of the main resonance O = w by using Melnikov 's 
method. Their insufficient analysis of the Melnikov function (in particular, that 
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of the role played by the initial phase 9) led them also into gross mistakes in their 
final conclusions. On the contrary, it was demonstrated by Chacon (19956) that 
the effect of the above parametric excitation in (1.2) for the general case of sub- 
harmonic resonances (i7 = puj, p an integer) is to suppress the chaotic behavior 
when a suitable initial phase is used and only for certain ranges of its ampli- 
tude. It was also shown (Chacon 19956) for the first time that such suitable 
initial phases are compatible with the surviving natural symmetry under the 
parametric excitation. It was also conjectured (Chacon 1998) that such max- 
imum survival of the symmetries of solutions from a broad and relevant class 
of systems, subjected both to primary chaos-inducing and chaos-suppressing 
excitations, corresponds to the optimal choice of the suppressory parameters; 
specifically, to particular values of the initial phase differences between the two 
types of excitations for which the amplitude range of the suppressory excita- 
tion is maximum. Rajasekar (1993) applied Melnikov's method to study the 
suppression of chaos in the Duflang-van der Pol oscillator 

X - c?x -\- ^x^ = - x^) X + f cos (wi) -I- ?7 cos (Ot -h O^) , (1.3) 

where the additional forcing r] cos (fit -|- O^) is the chaos-suppressing excitation, 

for the single case of the main resonance fl = o). He pointed out the relevant 
role of the initial phase (of the chaos-suppressing excitation) on the suppression 
scenario for the first time, and he also deduced the analytical expression of the 
boundaries of the regions in the rj — cj) phase plane where homoclinic chaos is 
inhibited. A generalization of Rajasekar 's approach concerning the relative ef- 
fectiveness of any two weak excitations in suppressing homoclinic/heteroclinic 
chaos is discussed in the work of Chacon (2002). There, general analytical ex- 
pressions are derived from the analysis of generic Melnikov functions providing 
the boundaries of the regions as well as the enclosed area in the amplitude-initial 
phase plane of the chaos-suppressing excitation where homoclinic/heteroclinic 
chaos is inhibited. Also, a criterion based on the aforementioned area was de- 
duced and shown to be useful in choosing the most suitable of the possible 
chaos-suppressing excitations. Cicogna & Fronzoni (1993) analyzed the Mel- 
nikov function associated with the family of systems 

X = f{x) - Sx + "f cos (iot) + eg (x) cos (fit + 9) , (1.4) 

where eg {x) cos (ilt -|- 9) is the chaos-suppressing excitation, for the single case 
of the main resonance O = w. They deduced both the suitable suppressory val- 
ues of the initial phase 9 and the associated chaotic threshold function (7/(5)^^ 
when the chaos-suppressing excitation acts on the system. General results 
(Chacon 1999) concerning suppression of homoclinic/heteroclinic chaos were 
derived on the basis of Melnikov's for the family (1.1) for the general case of 
subharmonic resonance (f2 = puj, p an integer). There, a generic analytical 
expression was deduced for the maximum width of the intervals of the initial 
phase Q for which homoclinic/heteroclinic bifurcations can be frustrated. It 
was also demonstrated that {0, 7r/2, tt, 37r/2} are, in general, the only optimal 
values of such initial phase, in the sense that they allow the widest ampli- 
tude ranges for the chaos-suppressing excitation. The work of Chacon (2001a) 
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presents general results concerning enhancement or maintenance of chaos for 
the family (1.1), where the connection with the results on chaos suppression 
was discussed in a general setting. It was also demonstrated that, in general, 
a second harmonic excitation can reliably play an enhancer or inhibitor role by 
solely adjusting its initial phase. The work of Chacon (20016) provides a pre- 
liminary Melnikov-method-based approach concerning suppression of chaos by 
a chaos-suppressing excitation which satisfies an ultrasubharmonic resonance 
condition with the chaos-inducing excitation. This approach was further ap- 
plied to the problem of the inhibition of chaotic escape from a potential well by 
incommensurate escape-suppressing excitations (Chacon & Martinez 2002). 



2 Basic theoretical approach 

To illustrate the theoretical approach with a paradigmatic example, consider a 
single Josephson junction subjected to a nonlinear dissipative term and driven 
by two harmonic excitations (Chacon et al. 2001) 

'x + smx = -a{l + ^cosx)x Fsin(wt) + (iFsm{nt + 9) , (2.1) 

where x and time are dimensionless variables, and x is proportional to the 
difference of potential between the two superconductors. It is also assumed that 
the terms of dissipation and excitation are regarded as weak perturbations and 
(iF sm{ilt + '^) is the chaos-suppressing excitation. The nonlinear dissipative 
term appears in the study of a single Josephson junction when the conditions 
are such that the interference effects between the pair and quasiparticle currents 
should be taken into account (Barone & Paterno 1982). The application of the 
Melnikov method to (2.1) yields the Melnikov function 

M± {to) = -C±A sin {uto) ± B sin (^o + B) , (2.2) 

with 



C = 8a (1 + 7/3), 
A = 27rFsech(^), 

B = 27r/3sech(^^^ . (2.3) 

Turning to the general case (1.1), let us assume that such a family of systems 
satisfies the requirements of the Melnikov method. Then, the application of the 
method to (1.1) provides the generic Melnikov function 

^th' (^o) = D± Ahar (wtq) + Bhar' (^to + , (2.4) 

where har (r) means indistinctly sin (r) or cos (r) , and A is a non-negative 
function, while D, B can be non-negative or negative functions, depending upon 
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the respective parameters for each specific system. In particular, D contains the 
effect of the damping, time-delay terms, and constant forces. In the absence of 
time-delay terms and constant drivings, D < 0, while one has the three cases 
D ^ when a constant driving and a time-delay term act on the system. Also, A 
and B contain the effect of the chaos-inducing and chaos-controlling excitations, 
respectively. Note that changing the sign of B is equivalent to having a fixed 
shift of the initial phase: B —B ^hh' ^ ^hh' ^ where the two 
signs before n apply to each of the sign superscripts of 4*. Therefore, B will be 
considered (for example) as a positive function in the following. As phase and 
initial time tq are not fixed, one may study the simple zeros of Mj^f^, (to) by 
choosing quite freely the trigonometric functions in (2.4). Therefore, consider, 
for instant, the Josephson junction given by (2.1). It is worth noting that 
the Melnikov functions (tq), (to) (cf. (2.4) and (2.2), respectively) 

are connected by linear relationships which are known for each specific system 
(1.1): 

To = To(to,w), (2.5) 

Therefore, the control theorems associated with any Melnikov function M^^, (tq) 
can be straightforwardly obtained from those associated with {to) (Chacon 
1999). 

2.1 Suppression of chaos 

As is well-known, the Melnikov method provides estimates in parameter space 

for the appearance of homoclinic (and heteroclinic) bifurcations, and hence for 
transient chaos. This means that in most of cases only necessary conditions for 
steady chaos (strange chaotic attractor) are obtained from the method. There- 
fore, one may always get sufficient conditions for the inhibition of even tran- 
sient chaos (frustration of homoclinic/heteroclinic bifurcation) and, a fortiori, 
for the inhibition of the steady chaos that ultimately arises from such a homo- 
clinic/heteroclinic bifurcation. This is the principal foundation of the utility of 
Melnikov method in predicting the suppression of (steady) chaos when a homo- 
clinic/heteroclinic bifurcation occurs prior to its emergence. For the Josephson 
junction (2.1) one has the following theorem (Chacon et al. 2001): 

Theorem 1 Let ^} = puj, p an integer, such that, for M+ (to) (M~ {to)), p = 

'^"^^in+i^^ {p ^ satisfied for some integers m and n. Then 

M"^ {to) always has the same sign, specifically M"^ {to) < 0, if and only if the 
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following condition is satisfied: 



/^min ^ P ^ /^max' 




max 2 ' 



cosh {Trn/2) 
cosh {ttcl)/2) ' 



(2.6) 



Now, the fohowing remarks are in order. 

First, one can test the suppression theorem theoretieahy by considering the 
limiting Hamiltonian case (a = 0). Notice that, in the absence of dissipation, 
(2.6) must be rewritten as ^ /3 ^ (3^^^, (3^-,^ = R, f3,,,^^ = R/p"^, since 

cannot now be zero. Thus, one obtains (Chacon et al. 2001) il = uj, Q = tt, 
and /? = 1 as a necessary and sufficient condition for suppressing stochasticity. 
(This resuh can be trivially obtained, to first perturbative order, from (2.2), 
(2.3) with a = 0, i.e., having M± (to) = for ah to- 
Second, the lower threshold for the chaos-suppressing amplitude, P^i^, takes 
into account the strength of the initial chaotic state through the factor (1 — C/A), 
since one usually finds that the corresponding maximal Lyapunov exponent 
increases as the ratio C/A decreases over a certain range of parameters. There- 
fore, for fixed-chaos inducing and chaos-suppressing frequencies (and hence R 
fixed), one would expect that will increase as A"*" is increased. Note that 
the corresponding upper threshold, /3max> does not verify this important prop- 
erty, which is because fH^^^^ arises from a necessary condition for the necessary 
condition yielding f3^^^ to be also a sufficient condition. This means that /3^ax 
slightly underestimates the upper threshold for the chaos-suppressing ampli- 
tude, as is numerically and experimentally observed in difi^erent instances. It is 
worth noting that this remark holds for any Melnikov function (2.4). 

Third, the asymmetry between the upper and lower homoclinic orbits (cf. 
(2.2), (2.3)) gives rise to two distinct sets of optimal initial phases that are 
suitable for suppressing chaos. The optimal suppressory values of 6 (hereafter 
denoted as Qopt) are those values allowing the widest amplitude ranges for the 
chaos-suppressing excitation (the use of the adjective is justified below in the 
discussion of the suitable intervals of initial phase difference for taming chaos) . 
Indeed, Theorem 1 requires having O = Qopt ^ tt, 7r/2, 0, 37r/2 (tt, 37r/2, 0, 7r/2) 
for p = An — 3,An — 2,An — 1, 4n (n = 1, 2, ...), respectively, in order to inhibit 
chaos when one considers orbits initiated near the upper (lower) homoclinic 
orbit. These distinct values are those compatible with the surviving natural 
symmetry under the additional forcing. Indeed, the dissipative, harmonically 
driven Josephson junction (/3 = 0) is invariant under the transformation 



X 



t 



t + 



(2n + 1) TT 



(2.7) 
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where n is an integer, i.e., if is a solution of (2.1) with /3 = 0, then 

so is [—x{t + (2n + 1) tt/w), —x{t + (2n + 1) tt/w)] . This pair of solutions may 
be essentially the same in the sense that they may differ by an integer number 
of complete cycles, i.e., 

x{t) = -x[t+{2n+l)Tr/uj]+2nl, (2.8) 

with I an integer, and they are termed symmetric. Otherwise, the time-shifting 
and sign reversal procedure yields a different solution, and the two solutions are 
termed broken-symmetric. When /3 > and is arbitrary the aforementioned 
natural symmetry is generally broken. The reason for that breaking is 

sin {m + @)y^ sin [m + e + (2n + 1) 7rfi/w] , (2.9) 

for arbitrary w, 17, and 8. Assuming a resonance condition O = poj, the survival 
of the above symmetry implies 

sin {pLot + &) = (-1)^+^ sin {puit + Q) . (2.10) 

Obviously, this is only the case for p an odd integer. For p an even inte- 
ger, one has the new transformation [x —x, t ^ t + {2n + 1)^/0), Q —^ 
© ± tt]. In other words, if [a;(f), a;(t)] is a solution for a value ©, then so is 
[-x{t + (2n + 1) tt/lu), -x{t + (2n 4- 1) tt/w)] for 9 ± tt. Thus, this explains 
the origin of the differences between the corresponding (at the same resonance 
order) allowed Qopt values for the two homoclinic orbits. Similar results have 
been found for the damped, driven pendulum mounted on a vertically oscillating 
point of suspension (Chacon 19956). Therefore, this maximum symmetry prin- 
ciple appears to be the common background in the mechanism of regularization 
by the application of resonant excitations. 

Fourth, the width of the allowed interval ]/3niin, /3max] fo'^ regularization is 



with R given by (2.6). Since i? is a positive function, there always exists a 
■maximum resonance order Pmax for suppression of homoclinic (and heteroclinic) 
chaos, for each fixed initial chaotic state (i.e., C/A fixed), which is 



-1/2- 



(2.12) 



where the brackets indicate integer part. From (2.12), one sees that Pmax in- 
creases as the ratio C/A is increased, which would associated with the decrease 
of the corresponding maximal Lyapunov exponent over a certain range of pa- 
rameters. For a given set of parameters satisfying the above theorem's hy- 
pothesis, as the resonance order p is increased the allowed interval ]/3inin> /^max] 
shrinks quickly for low frequencies. This means that initial chaotic states cannot 
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necessarily be regularized to periodic attractors of arbitrary long period, since 
numerical experiments indicate that the regularized responses are typically a 
period-1 attractor for p = 1 and a period-2 attractor for p = 2. On the other 
hand, the asymptotic behavior A/3 (lu oo) = oo (the remaining parameters 
being held constant) means that chaotic motion is not possible in this limit, as 
expected. 

Fifth, to establish the suppression theorem corresponding to any Mclnikov 
function (2.4), it is enough to transform M^^, (tq) into the form given by (2.2). 
Therefore, taking into account (2.5) and the aforementioned Qopt values, one 
finds that in general there exist at most four suitable optimal values for the 
suppressory initial phase difference between the two (commensurate: fl = poj) 
excitations: 0, 7r/2, tt, 37r/2. 

Sixth; It has been stated above that the suitable values of the initial phase 
difference (between the two excitations involved) given by Theorem 1 are op- 
timal, in the sense that they allow the widest amplitude ranges for the chaos- 
suppressing excitation. One therefore could expect reliable control of the dy- 
namics over certain suitable phase difference intervals, which would be centered 
on such optimal values, although this would imply a reduction of the respective 
amplitude ranges. It has been deduced (Chacon 1999; Chacon et al. 2001) that 
there always exists a maximum-range interval 

[Qopt ^ AOmaxj Qopt + ABnjax] i (2-13) 

of permitted initial phase differences for homo clinic/ hetero clinic chaos inhibi- 
tion, where 

AQmax = arcsin ( -r ) . (2.14) 



For each value of Q belonging to this interval there exists a reduced interval 
(with regard to the limiting case where the only suitable values of G are Qopt) 
of amplitudes of the chaos-suppressing excitation which is 

/3min (Q = ©opt ± Ae) < /3 (6 = e„j,t ± Ae) , 
/?min (0 = Qopt ± Ae) = (^i-^^i?sec(Ae), 

/3:„ax (Q = Qopt ± Ae) = (2.15) 

where R is given by (2.6) and ^ AG ^ AQmax- Thus, the width of the range 
for the chaos-suppressing amplitude is 

i.e., for fixed C, A and AQ there always exists a maximum resonance order Pmax 
for homo clinic chaos suppression which is 



Pmax (Q = Qopt ± Ae) 



cos (AQ) 



-y/l - sin(Aeinax) 



(2.17) 
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where the brackets indicate integer part. Also, one can put 



AG, 



5-° 



(2.18) 



Thus, one can use a linear approximation for A6max suitable for chaotic motions 
arising away from the limiting case of tangcncy between the stable and unstable 
manifolds (C/^ <C 1). It is worth mentioning that the last inequality is usually 
associated with the observation of steady chaos (strange chaotic attractor). 



2.2 Enhancement of chaos 

It has been mentioned above that the mechanism for suppressing homoclinic 
(and heteroclinic) chaos is the frustration of a homoclinic/heteroclinic bifurca- 
tion, which prevents the appearance of horseshoes in the dynamics. Chacon 
(2001a) showed that the enhancement of the initial chaos is achieved by moving 
the system from the homoclinic tangency condition even more than in the initial 
situation with no second periodic excitation. Thus, enhancement of chaos can 
mean increasing the duration of a chaotic transient, passing from transient to 
steady chaos, or increasing the maximal Lyapunov exponent. Consider again 
that the family of systems modeled by (1.1) satisfies the requirements of the 
Melnikov method. Similarly to the preceding discussion of the suppression of 
chaos, one can assume any particular form of the Melnikov function (2.4) to dis- 
cuss the enhancement of chaos. Therefore, consider, for instance, the following 
nonlinearly damped, biharmonically driven, two-well Duffing oscillator: 

x-x + (3x^ = -dx -h Fcos {ut) - rjpx^ cos {fit + 6) , (2.19) 

where ri,fl, and 8 are the normalized amplitude factor, frequency, and initial 
phase, respectively, of the chaos-controlling parametric excitation (0 < 77 ^ 1), 
and /3, S, n, F, and to are the normalized parameters of the potential coefficient, 
damping coefficient, damping exponent, chaos-inducing amplitude, and chaos- 
inducing frequency, respectively (0 < 6 , F 1, (3 > 0, n = 1,2, ...). The appli- 
cation of the Melnikov method to (2.19) yields the Melnikov function 

M=^(io) = D±Asm (wio) - C^sin (Oi + 6) , (2.20) 

with 

A ^ 0) .c.Fsech(^), 

C ^ l^{Q'^ + 4i}^)csch^^y (2.21) 
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where the positive (negative) sign refers to the right (left) homochnic orbit of the 
underlying integrable two- well DufSng oscillator {S = F = r] = 0), and where 
B (m, n) is the Euler beta function. It has been demonstrated (Chacon 2001a) 
that, in general, a second harmonic excitation can reliably play an enhancer 
or inhibitor role solely from adjusting its initial phase. The Melnikov function 
M+(to) will be used here to illustrate the approach to the enhancement of 
chaos. Indeed, consider that, in the absence of any second parametric excitation 
(C = 0), the associated Melnikov function Mg'"(fo) = — l^l -|->lsin (cuto) changes 
sign at some to, i.e., < ^. If one now lets the second excitation act on 
the system such that C ^ ^ — jDj, this relationship represents a sufficient 
condition for M+(<o) to change sign at some to- Thus, a necessary condition 
for M+{to) to always have the same sign {M+{to) < 0) is C > A - |Z)| = 
Cmin- It was above mentioned (Chacon 1999) that a sufficient condition for 
C > Cmin to also bc & Sufficient condition for inhibiting chaos is Q = pu) 
(subharmonic resonance condition), C < Cmax = ^./p^, p an integer, and that 
-^o'(io) and — Cmin. max sin (fito + ©) are in opposition. This condition yields 
the optimal suppressory values = ©opt- It was demonstrated (Chacon 

2001a) that imposing MQ{to) to be in phase with —Cmin, max sin {^to + O) is a 
sufficient condition for M~^{to) change sign at some to- This condition provides 
the optimal enhancer values of the initial phase, Qlpt, in the sense that A'/+(to) 
presents its highest maximum at Qlpt^ i-^-j obtains the maximum gap from 
the homoclinic tangency condition. Now, the following remarks are in order. 

First, for a given homoclinic orbit forming (part of) a separatrix, one has in 
general (i.e., for any Melnikov function (2.4)) that 

K7t-®Tpt\=^, (2-22) 

for each resonance order. 

Second, for C = Cmin there always exists a maximum-range interval 

[QTpt - Ae-t (C = Cmin) , + AG-^, (C = Cmin)] (2.23) 

of permitted initial phases for enhancement of chaos in the sense that, for values 
of © belonging to that interval, the maxima of M+(to) are higher than those of 
M^(to)- Similarly, for C = Cmax there always exists a different maximum-range 
interval 

[ei;^ - A0St (C = Cmax) , Qi;^ + A©SLx (C = Cma.)] (2-24) 

of allowed initial phases for enhancement of chaos, and also 

A0-^ (C = Cmax) > A©-^, (C = Cmin) , (2-25) 

which is a consequence of the dissipation. It must be emphasized that the 
definition of Q^pt is general; i.e., it refers to any resonance and any Melnikov 
function (2.4). 

Third, for general scparatrices, i.e., those formed by several homoclinic and 
(or) heteroclinic loops, the above scenario of control of chaos holds for each 
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homoclinic (hctcroclinic) orbit. However, it is common to find that the dif- 
ferent homoclinic (heterochnic) orbits of a given separatrix yield distinct Qlpf 
values. This is a consequence of the survival of the symmetries existing in the 
absence of the second excitation. Thus, the actual scenario is usually more 
complicated. For instance, let O^yj^^, 0™^; be the optimal values associated 
with the right and left homoclinic orbits, respectively, of a typical separa- 
trix with a "figure-of-eight" loop, as in the two-well Duffing oscillator (2.19). 
One then obtains that the best chance for enhancing chaos should now be at 

©opt ~ {^opt,r - ©opm) /2 mod (27r) . See Chacon (2001a) for more details. 
2.3 Further developments 

The case of subharmonic resonance between the chaos-inducing and chaos- 
controlling frequencies has been briefly discussed above. However, a number of 
theoretical (Salerno 1991; Salerno & Samuelsen 1994), numerical (Braiman & 
Goldhirsch 1991), and experimental (Uchida et al. 1998) studies show that chaos 
can be reliably controlled by other non-subharmonic resonances. The work of 
Chacon (20016) presents a Melnikov-method-based approach concerning reduc- 
tion of homoclinic and heteroclinic instabilities for the family of systems (1.1) 
where the harmonic excitations verify an ultrasubharmonic resonance condition: 
i}/ui = p/q, q > 1 {p q), p,q positive integers and O (w) the chaos-suppressing 
(inducing) frequency. Such general results can be used to approach the case 
of incommensurate chaos-suppressing excitations by means of a series of ever 
better rational approximations, which are the successive convergents of the infi- 
nite continued fraction associated with the irrational ratio Q/uj. This procedure 
has been much employed in characterizing strange non-chaotic attractors in 
quasipcriodically forced systems as well as in studying phase-locking phenom- 
ena in both Hamiltonian and dissipative systems. To illustrate the method one 
intentionally chooses the golden section Q/oj = <I> = (-\/5 — l) /2, since it is the 
irrational number which is the worst approximated by rational numbers (in the 
sense of the size of the denominator) . As is well-known, the golden section can 
be approximated by the sequence of rational numbers (f^/w), = Fi-i/Fi where 
Fi = 1, 1, 2, 3, 5, 8, are the Fibonacci numbers such that limj^oo = ^• 

For each {Q/u))^ one replaces each quasiperiodically excited system 

dU(x) 

X -\ = -d(x,x) + pJx,x)har{ojt) +Ps (x,x) har' (^ojt + h') (2.26) 

dx \ y ' 

by the respective periodically excited system 
dU{x 



x + 



= -d{x,x) + Pc{x,x)har{ojt) +Ps {x,x) har' ^"^T"^* + 



dx 

(2.27) 

giving a sequence of periodically excited systems whose associated frequen- 
cies satisfy an ultrasubharmonic resonance condition. The work of Chacon & 
Martinez (2002) applied this approach to the problem of the reduction of chaotic 
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escape from a potential well using the simple model 

x = x- I3x^ -Sx + ^ sin (ujt) - Prjx'^ sin (Of + 6) , (2.28) 

where (3r]x'^ sin {fit + 6) is the escape-suppressing excitation. They found that, 
for irrational escape-suppressing frequencies, the effective escape-reducing initial 
phases are found to lie close to the accumulation points of the set of suitable 
initial phases that are associated with the complete series of convergcnts up to 
the convergent giving the chosen rational approximation. 

A Melnikov-method-based approach (Chacon 2002) was presented concern- 
ing the relative effectiveness of harmonic excitations in suppressing homoclinic 
(and heteroclinic) chaos of the family (1.1) for the main resonance between 
the chaos-inducing and chaos-suppressing excitations. A criterion based on the 
area in the supprcssory amplitude/initial phase parameter plane, where suppres- 
sion of homoclinic chaos is guaranteed, was deduced and shown to be useful in 
choosing the most suitable of the possible chaos-suppressing excitations. Addi- 
tionally, the choice of the most suitable chaos-suppressing excitation was shown 
to exhibit sensitivity to the particular initial chaotic state. 

The work of Chacon et al. (2003) presents general findings concerning control 
of chaos for the family 

— = -d {x, x) hch,i {x, x) i^c/i,i(i)+X] ^''"'i ^) -^'^o-jW' (2-29) 

where U (x) is a general potential, —d (x, x) represents a generic dissipative 
force, X^^j hch,i Fch.i{t) is a general multiple chaos-inducing excitation, 

and J2j=i h J {t) is an as yet undetermined suitable multiple chaos- 

controlling excitation, with Fch,i{t), Fco,j{t) being harmonic functions of com- 
mon frequency uj and initial phases (i = 1,...,N), ipj (j = 1,...,M). The 
effectiveness of this approach in suppressing spatio-temporal chaos of chains 
of identical chaotic coupled oscillators was demonstrated through the example 
of coupled Duffing oscillators, where coherent oscillations were achieved under 
localized control. 

The work of Chacon et al. (2002) studied the robustness of the suppression 

of bidirectional chaotic escape of a harmonically driven oscillator from a quartic 
potential well by the application of weak parametric excitations. It was numer- 
ically shown that Melnikov-method-based theoretical predictions also work for 
harmonic escape-inducing excitations in the presence of external noise, and for 
chaotic-escape-inducing excitations having a sharp Fourier component with a 
sufficiently high power. 

The method proposed in the work of Lenci & Rega (2003) consists of choosing 
the shape of external and/or parametric periodic excitations, which permits 
one avoid, in an optimal manner, a homoclinic bifurcation. They numerically 
investigated the effcctivciness of the control method with respect to the basin 
erosion and escape phenomena of a perturbed Helmholtz oscillator. 
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3 Some applications 



3.1 Taming chaotic escape from a potential well 

The work of Chacon et al. (1996, 1997, 2001) and Balibrea et al. (1998) 

applies the above Melnikov-method-based approach to the problem of chaotic 
escape from a potential well. This is a general and ubiquitous phenomenon in 
science. Indeed, one finds it in very distinct contexts: the capsizing of a boat 
subjected to trains of regular waves (Thompson 1989), the stochastic escape of 
a trapped ion induced by a resonant laser field (Chacon & Cirac, 1995), and 
the escape of stars from a stellar system (Contopoulos et al. 1993) are some 
important examples. Remarkably, such complex escape phenomena can often 
be well described by a low-dimensional system of difierential equations. The 
case considered by Chacon and coworkers is that where escape is induced by an 
external periodic excitation added to the model system, so that, before escape, 
chaotic transients of unpredictable duration (due to the fractal character of the 
basin boundary) are usually observed for orbits starting from chaotic generic 
phase space regions (such as those surrounding separatrices) , in both dissipative 
and Hamiltonian systems. In particular, Chacon et al. (1996) studied the 
simplest model for a universal chaotic escape situation: 

.-..,.^=-....s.M,.(-;-:-<^-«)), (3.) 

where /3rya;^ sin (fit + 6) and /?rya;sin (Ot + 6) are the (independently consid- 
ered) escape-suppressing parametric excitations. It was demonstrated that the 
parametric excitation of the linear (quadratic) term suppress chaotic escape 
more efficiently than that of the quadratic (linear) term for small (large) driv- 
ing periods of the primary chaos- inducing excitation. Chacon et al. (1997) 
studied the inhibition of chaotic escape of a driven oscillator from the cubic 
potential well that typically models a metastable system close to a fold: 

'x + x-l3x^ = -( 2^^^ \ +^cos {ujt) - 'qx COS {Vlt + Q) , (3.2) 

where 8ix and {52x'^ + dzx'^^ x are the (independently considered) linear and 
nonlinear damping terms, respectively. They demonstrated that the effective- 
ness of a parametric excitation at suppressing chaotic escape from such a cubic 
potential well diminishes as the system approaches a period- 1 parametric reso- 
nance^ and that, for linear damping, the parametric excitation inhibits chaotic 
escape more efficiently than for nonlinear damping. The role of a nonlinear 
damping term, proportional to the nth power of the velocity, on the escape- 
inhibition scenario is considered in the work of Chacon et al. (2001): 

X + X — x^ = —Sx\x^ ^ +^cos {ujt) +r]x^ cos {Q.t + Q) , (3.3) 

where r]x'^ cos {Q,t + 9) is the escape-suppressing parametric excitation. In this 
case, the effectiveness of the parametric excitation of the quadratic potential 
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well at inhibiting chaotic escape diminishes as the system approaches either 
a period-1 or a period-2 parametric resonance. Also, the effectiveness of the 
parametric excitation in the presence of the nonlinear dissipative force is less 
than for a linear dissipative force. 

3.2 Taming chaotic solitons in Prenkel-Kontorova chains 

Control of chaos in spatially extended systems is one of the most important 
and challenging problems in the field of nonlinear dynamics. Instances of possi- 
ble applications include the stabilization of superconducting .Josephson-junction 
arrays (Barone & Paterno 1982), periodic patterns in optical turbulence, and 
semiconductor laser arrays (SchoU 2001), to cite just a few. Martinez & Chacon 
(2004) presented a Melnikov-method-based general theoretical approach to con- 
trol chaotic solitons in damped, noisy and driven Frenkel-Kontorova chains. 
Specifically, they studied the model 

K 

Xj + — sin {2nXj) = a^j+i — 2xj + xj^i — axj + F cos {ivt) 

+PFcos{nt + ^)+^{t), (3.4) 

where l3F cos (Ot + (p) is the chaos-suppressing excitation, and ^ (t) is a bounded 
noise term. They obtained an effective equation of motion governing the dynam- 
ics of the soliton center of mass for which they deduced Melnikov-method-based 
predictions concerning the regions in the control parameter space where homo- 
clinic bifurcations are frustrated. Numerical simulations indicated that such 
theoretical predictions can be reliably applied to the original Prenkel-Kontorova 
chains, even for the case of localized application of the soliton-taming excitations. 
It is worth mentioning that the same effectiveness of such a localized control 
in suppressing spatio-temporal chaos of chains of identical chaotic coupled os- 
cillators was demonstrated through the example of coupled DufBng oscillators 
(Chacon et al. 2003). 

3.3 Taming chaotic chcirged particles in the field of an 
electrostatic wave packet 

The interaction of charged particles with an electrostatic wave packet is a basic 
and challenging problem appearing in diverse fundamental fields such as astro- 
physics, plasma physics, and condensed matter physics. While the Hamiltonian 
approach to this problem is suitable in many physical contexts, the considera- 
tion of dissipative forces seems appropriate in diverse phenomena such as the 
stochastic heating in the dynamics of charged particles interacting with plasma 
oscillations. In any case, stochastic (chaotic) dynamics already appears (can 
appear) when the wave packet solely consists of two electrostatic plane waves. 
Such a non-regular behavior of the charged particles may yield undesirable ef- 
fects on a number of technological devices such as the destruction of magnetic 
surfaces in tokamaks. Thus, apart from its general intrinsic interest, the problem 
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of rcgularization of the dissipativc dynamics of charged particles in an electro- 
static wave packet by a small-amplitude uncorrelated wave (which is added to 
the initial wave packet) is especially relevant in plasma physics. Chacon (2004) 
considered the simplest model equation to examine this problem: 

6 

it + 6x = [Eq sin (kox — cjot) + Ec sin (kcX — tOct)] 

m 

Es sin (kgX — Wst) , (3.5) 

m 

where E^ sin {kcX — Uct) and Eg sin [kgX — uj^t) are the chaos- inducing and chaos- 
suppressing waves, respectively. In a reference frame moving along the main 
wave E{) sin {k^x — loqI) , (3.5) transforms into a perturbed pendulum equation 
which is capable of being studied by means of Melnikov's method. Two sup- 
pressory mechanisms were identified: One mechanism requires chaos-inducing 
and chaos-suppressing waves to have both commensurate wavelengths and com- 
mensurate relative phase velocities, while the other allows chaos to be tamed 
when these quantities are incommensurate. 



4 Conclusions and open problems 

The present review summarizes some of the main results and applications of a 
preliminary theoretical approach to control chaos in dissipativc, non- autonomous 
dynamical systems, capable of being studied by Melnikov's method, by means of 
periodic excitations. Diverse extensions and applications of the current theory 
remain to be developed. Among them: 

(i) To obtain the boundaries of the rcgularization regions in the control 
parameter space for the case of a general resonance (not just the main) between 
the involved excitations. 

(ii) To extend the theoretical approach to (some family of) multidimensional 
systems capable of being studied by (some generalized version of) Melnikov's 
method. 

(iii) To develop a multiharmonic control theory beyond the main resonance 
case. 

(iv) To extend the theoretical approach for the case of periodic excitations 
to the case of random excitations. 

(v) To obtain analytical approximations of the regularized responses for the 
deterministic case of a general resonance between the chaos-inducing and chaos- 
suppressing excitations. 

(vi) To extend the current theory described for harmonic excitations to the 
case of general periodic excitations (both chaos- inducing and chaos-controlling). 
In particular, the waveform effect should be taken into account in the control 
scenario. 

(vii) To extend the current theory to the case where the chaos-controlling 
excitation is a parametric excitation of the amplitude of the chaos-inducing 
excitation, as well as to the case where it is a parametric excitation of the 
frequency of the chaos-inducing excitation. 
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(viii) To apply the current theory to ratchet systems to improve the directed 
energy transport. 

(ix) To apply the current theory to control chaotic population oscillations be- 
tween two coupled Bosc-Einstein condensates with time-dependent asymmetric 
potential and damping. 

Acknowledgement 2 The author acknowledges financial support from. Span- 
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